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Abstract: We employ the lattice QCD data on Taylor expansion coefficients to extend our previous paramctrization 

of the equation of state to finite baryon density. When we take into account lattice spacing and quark 
mass dependence of the hadron masses, the Taylor coefficients at low temperature are equal to those of 
hadron resonance gas. Parametrized lattice equation of state can thus be smoothly connected to the hadron 
resonance gas equation of state at low temperatures. 
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One of the methods to extend the lattice QCD calculations to non-zero chemical potential is Taylor expansion. 
In that approach pressure is Taylor expanded in chemical potentials, and the Taylor coefficients are calculated on 
the lattice at zero chemical potential. In this contribution we use the results of the most comprehensive lattice 
QCD analysis of the Taylor coefficients to date [1, 2] to construct a parametrization of an equation of state (EoS) 
for finite baryon density. As in our earlier parametrization of the EoS at zero chemical potential [3], we require 
that our parametrization matches smoothly to the hadron resonance gas (HRG) at low temperatures. 
Taylor coefficients are derivatives of pressure with respect to baryon and strangeness chemical potential: 

1 rpi+j Qi Qj 

c.(T) = ^ — ^ — P(T,, fl =0,, s =0). (1) 
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As we discussed in [4], all the coefficients evaluated in Refs. [1, 2] are well below the HRG results. This discrepancy 
can largely be explained by the lattice discretization effects on hadron masses: When the hadron mass spectrum 
is modified accordingly (for details see [5]), the HRG model reproduces the lattice data, see Fig. 1 of Ref. [4]. 
Interestingly this change can be accounted for by shifting the modified HRG result of purely baryonic coefficients 
towards lower temperature by 30 MeV. The situation is similar for other Taylor expansion coefficients [5] , although 
the strange coefficients might favour slightly smaller shift. Based on this finding and because the lattice data agree 
so well with the modified HRG we suggest that cutoff effects can be accounted for by shifting the lattice data by 
30 MeV. We show the effect of such a shift in the left panel of Fig. 1, where we plot the HRG curve with physical 
masses (dashed line) and compare it with the lattice data, where all the points below 206 MeV temperature are 
shifted by 30 MeV, and the 209 MeV point by 15 MeV. For further confirmation of this procedure we also plot 
the recent HISQ result of C20 [6] in Fig. 1 (right): At low temperatures the shifted p4 data agree with the HISQ 
data. 

We parametrize the shifted data using an inverse polynomial of four (C20), five (en and C02), or six (fourth and 
sixth order coefficients) terms: 

m 

c ^) = £^g-+4 B I (2) 

k — l 

where cf^ is the Stefan-Boltzmann value of the particular coefficient, and the powers n^y are required to be 
integers between 1 and 42. T = (T — T s )/R with scaling factors T 3 — and R = 0.15 GeV for the second order 
coefficients and T 3 =0.1 GeV and R = 0.05 GeV for all the other coefficients. We match this parametrization 
to the HRG value at temperature T$w = 155 MeV by requiring that the Taylor coefficient and its first, second, 
and third derivatives are continuous. Since the recent lattice data obtained using HISQ action [6] shows that the 
second order coefficients approach their Stefan-Boltzmann limits slowly, we require that their value is 95% of their 
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Fi gUf6 1. The parametrization (solid line) and HRG value (dashed) of the C20 (left) and C40 (right) coefficients compared with 
the shifted p4 data (see the text). The recent lattice result for C20 with the HISQ action [6] is also shown. The arrows 
depict the Stefan-Boltzmann values. 
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Figure 2. (Left) The parametrization (solid line) and HRG value (dashed) of the cqq coefficient compared with the shifted p4 
data (see the text). 

(Right) Speed of sound as function of temperature on various iscntropic curves with constant entropy per baryon. 



Stefan-Boltzmann value at 800 MeV temperature. These constraints fix four (or five) of the parameters atij . The 
remaining parameters, are fixed by a x 2 fit to the lattice data. As an example we show the parametrized C20, C40, 
and C60 in Figs, f and 2. 

Once the coefficients are known, pressure can be written as 

^ = E^^)(f) i (f) i ' ( 3 ) 

and all the other thermodynamical quantities can be obtained from Eq.(3) by using the laws of thermodynamics, 
his kind of an expansion breaks down at large fi. However, at baryon densities of interest here, the contribution 
from coefficients of particular order is clearly below (< 20%) the lower order contribution. Thus the approximation 
of the EoS in terms of fourth- and sixth order expansion looks reasonable. As pressure at /is = 0, i.e. the coefficient 
coo, we use our earlier parametrization s95p-vl [3]. We describe the EoS in the right panel of Fig. 2 by showing 
the speed of sound on various isentropic curves with constant entropy per baryon [7]. The curves at s/ub = 400, 
65, and 40 are relevant at collision energies ^/snn = 200, 39 and 17 GeV, respectively. At s/ub = 400 (dotted 
line), the EoS is basically identical to the EoS at hb — (thin solid line). At larger baryon densities the effect of 
finite baryon density is no longer negligible. The larger the density, the stiffer the EoS above, and softer below 
the transition temperature. 

Furthermore, additional structure begins to appear around the transition temperature with increasing density. 
We expect this structure to be an artifact of our fitting procedure: Our fit is too good and elevates errors to 
features leading to additional ripples in the speed of sound. Another interesting feature in the equation of state 
is the rapid change of the speed of sound around T sw = 155 MeV and another change around T ~ f 85 MeV. The 
latter has its origin in our baseline [is = EoS. It follows hadron resonance gas up to T — 184 MeV temperature 
causing a change in the speed of sound at that temperature. Nevertheless, when pressure is plotted as a function 
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of energy density, these structures are hardly visible. Therefore we do not expect them to affect the buildup of 
flow and the evolution of the system, and consider our parametrization a reasonable first attempt. 
We have studied the effect of the EoS on flow by calculating elliptic flow in Pb+Pb collision at the full SPS 
collision energy (^/snn = 17 GeV). Our results are similar to those we have shown earlier [4]: Even if in an ideal 
fluid calculation at RHIC energy proton V2(pr) is sensitive to the order of phase transition [8], at SPS energy 
both proton and pion V2(pr) are insensitive to it. 

To summarise, we have shown that a temperature shift of 30 MeV is a good approximation of the discretization 
effects in the lattice QCD data obtained using p4 action. We have constructed an equation of state for finite 
baryon densities based on hadron resonance gas and lattice QCD data. At the full SPS energy (^/snn = 17 GeV) 
the pT-differential elliptic flow is almost insensitive to the equation of state. This is bad news for the experimental 
search of the critical point, since a change from a first order phase transition to a smooth crossover does not cause 
an observable change in the flow. 
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